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Abstract. Starting from the Helmholtz free energy we calcu- 
late analytically first- and second-order derivatives, as inter- 
nal energy and specific heats, for the ideal system and the ex- 
change and correlation interactions covering a broad range of 
degeneracy and relativity. The complex physics of Coulomb 
interactions is expressed by Pade Approximants, which reflect 
the actual state of our knowledge with high accuracy. We as- 
sume complete ionization and provide a base system of ther- 
modynamical functions from which any other thermodynami- 
cal quantities can be calculated. We chose for the base system 
the free energy, the pressure, the internal energy, the isothermal 
compressibility (or density exponent), the coefficient of strain 
(or temperature exponent), and the isochoric specific heat. By 
means of the latter potentials entropy, isobaric specific heat 
and adiabatic temperature gradient can be determined. We give 
comparisons with quantities which are composed by numeri- 
cal second-order derivatives of the free energy and show that 
numerical derivatives of the free energy as calculated, for in- 
stance, from EOS tables, may produce discontinuities for astro- 
physically relevant quantities as, e.g., the adiabatic temperature 
gradient. Adiabatic temperature gradients are shown for differ- 
ent chemical compositions (hydrogen, helium, carbon). Finally 
the used formalism of Pade Approximants allows immediate 
incorporation of recent results from many particle statistics. 

Key words: Equation of state - Plasmas - Stars: interiors - 
Stars: evolution - Stars: low-mass, brown dwarfs - Stars: white 
dwarfs 



1. Introduction 

For an accurate modelling of stellar objects we have to compute 
a complete set of thermodynamical quantities which meets the 
physical conditions of various evolutionary stages. Based on 
the framework presented in Stolzmann & Blocker (1996a, here- 
after paper I) we derive further thermodynamical potentials in 



order to provide a base system from which any other thermo- 
dynamical quantity can be calculated. Such a system requires 
six quantities consisting of the Helmholtz free energy, two of 
its single derivatives, one of its second mixed derivatives, and 
two of its second pure derivatives with respect to temperature 
and density (see e.g. Dappen et al. 1988). We complement the 
Helmholtz free energy with the pressure and the internal en- 
ergy, with the isochoric specific heat, and with the temperature 
and density exponents in the equation of state (EOS), xt and 
Xp- This paper will be devoted to the ideal, exchange and coiTe- 
lation contributions according to the expansions given in paper 
I, viz. for case a: weak relativity and arbitrary degeneracy, and 
case b: strong degeneracy and arbitrary relativity. These cases 
cover a large area in the density-temperature plane. The cor- 
relations between the charged particles are formulated by the 
technique of Pade Approximants as in paper I. Here we give 
improved versions of our earlier applied Pade Approximants. 
In particular we have rearranged terms of the quantum virial 
function. 

Explicit expressions for the Helmholtz free energy and the 
pressure have already been derived in paper I for ideality 
as well as for exchange and coiTelation interaction. This pa- 
per aims at providing analytical formulae with high accuracy 
which supply quick computing with reliable accuracy in prac- 
tical applications. Furthermore the expressions of the thermo- 
dynamical functions given here can be easily included as a part 
of the thermodynamical description of arbitrary degree of ion- 
ization. 

This paper is organized as follows. In Sect. || we list the set 
of thermodynamical potentials to be considered. Sect. || gives 
a brief overview of the concept for the calculation of the EOS 
terms. Sect. ^ deals with the detailed determination of the ideal, 
exchange and coiTelation parts for the EOS. Numerical results 
and comparisons are presented in Sect. |], and a summary is 
given in Sect. |[ 
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2. Thermodynamical relations and identities 

We summarize briefly some well-known standard relations, 
which are frequently used to provide the thermodynamics for 
astrophysical applications. 

We have to determine the first and second-order quantities 
which are related by the ratio Cp/Cy of isochoric to isobaric 
specific heat: 



Cp V xl 
7 = — = IH — , 

Cv Kt^ CyT 



(1) 



The (isobaric) thermal expansion coefficient Ap can be ex- 
pressed by 



Kt 

Ap = — CvT7G = Kt^s , 

with 7g being the Griineisen coefficient 



PV 



-XT 



CyT' 

$s is flie coefficient of strain 

^"8 = PXT , 

and Kt is the isothermal compressibility 



(3) 



(4) 



(5) 



The so-called temperature and density exponents in the equa- 
tion of state (Cox & Giuli 1968) are defined by 



'd\nP\ _ T (dP\ 
^'^ ^ ydlnTJp ^ P [dfj 



V 



dlnP 
d\np 



P [dV 



(6) 



(7) 



The adiabatic temperature gradient defined by Vad = 
{d\nT/d\n P)s (S denotes the entropy) can be expressed by 



'ad 



PV , _ PV XT 
CpT ^ CpTxp 



(8) 



Another possibility to calculate the adiabatic gradient is given 
by the three adiabatic exponents Fi, and Fa (Cox & Giuli 
1968, Rogers et al. 1996) 



1 r, - 1 



'ad 



r2 Ti 

which can be obtained by 



Ti 



7 



7 



PKt 



7 - PKtJg 



(9) 



Ts = 1+7G .(10) 



In order to calculate the quantities given by Eqs. (ll|)-(y^ we 
start with the Helmholtz free energy F{V,T,N) and determine 
by means of standard thermodynamic relations the Gibbs en- 
ergy. 



(11) 



the pressure, 

PV = G-F , 

and the internal energy. 



U = F-T 



dT 



V,N 



(2) a 



leading to the entropy S according to 

TS - F 

and to the isochoric specific heat via 



dT 



{dry 



(12) 



(13) 



(14) 



(15) 



Moreover, we have to calculate the coefficient of strain and the 
inverse compressibility (bulk modulus) by means of Eqs. (^- 

(§■ 

3. Theoretical model 

As in paper I we start with the Helmholtz free energy of a 
fully ionized plasma consisting of ideal and Coulomb interac- 
tion parts 



F{T,V,Na)^}2^a'' + F' 



(16) 



where the Coulomb term contains the following parts: 

F^°^' = F^, + F^, + F^,+F.^ + F^, , (17) 

with X and c marking the exchange and correlation contribu- 
tions and cq the quantum correction term. The pairs ab denote 
the interaction between particles of species a and b (electrons, 
ions), respectively. 

For convenience we introduce dimensionless thermody- 
namical potentials defined by 



/ = 



F 



NkT ' 



9 = 



S 



1 



G 

NkT ' 

1 



U 



Nk ' kT Kt nkT 



P = 



nkT 



NkT 



(18) 



nkT' '-^^Wk^'^^ 



where n — N /V refers to the total particle number density of 
the ions or the electrons, and k denotes the Boltzmann con- 
stant. Eq. (|l6|), or correspondingly the potentials summarized 
by Eq. (|l8|) and Eq. can be written by 

S = N,kT[af+al^ + ai:\ + 

NkT [a^ + gI + a'r^ + a^] (20) 

where T. = {F, G, P ■ V, U, S ■ T, V/Kt, $s -V, Cy ■ T} 
and a = {f,g,p,u,s,l/kT,4>s,cy} symbolize the various 
thermodynamic functions defined in the previous section. 
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4. Thermodynamical potentials 

This section deals in detail with the calculation of the poten- 
tials listed in Eq. (^). By introducing definitions and for the 
sake of the integrity we repeat here few expressions concerning 
Helmholtz free energy and pressure, which are already given in 
paper I. 

4.1. Ideality 

The ideality of the nonrelativistic and nondegenerate ions 
is described by the well-known classical expressions for the 
Helmholtz free energy 



[ln(niAf)-l] 



(21) 



The thermal de Broglie wavelength for particles of species a is 
Aq = 2nfi / \/2ma'KkT . In the classical description we get for 
the pressure, the compressibilty, and the coefficient of strain 



1 

Lid 



= 1 



(22) 



and for the internal energy and the isochoric specific heat 
2 



(23) 



The ideal pressure of the electrons at any relativity and de- 
generacy (Chandrasekhar 1939, Cox & Giuli 1968) can be cal- 
culated by 



A3 



-^3/2(V',A) + -AJ5/2(V', A) 



(24) 



on the expressions given in Eqs. (|24|)-(p7[) by numerical inte- 
grations to determine the adiabatic temperature gradient. Re- 
cently, Blinnikov et al. (1996) and Miralles & Van Riper (1996) 
presented parametrizations and used various approximations to 
evaluate the fully relativistic ideality for the set of thermody- 
namical potentials, which are listed in Sect. ||. However, all 
calculational schemes are characterized by an immense effort 
in order to determine temperature- and density derivations even 
for asymptotic regions. For details, see e.g. Miralles & Van 
Riper (1996). 

We pursue to include relativistic effects over a broad re- 
gion of astrophysically relevant densities and temperatures for 
ideality and exchange. Furthermore we evaluate the set of ther- 
modynamical potentials analytically avoiding the well-known 
noise problems of second-order quantities apparent in purely 
numerical approaches. 

This can be realized by introducing two approximations: 

a) arbitrary degeneracy, but weak relativity {case a) and 

b) arbitrary relativity, but strong degeneracy {case b) 

A numerical study on the density-temperature validity region 
of these approximations was carried out in paper I. Note, that 
case a, based on A-expansions ^ 0{\'^), is limited to T ^ 
2 • lO^K and p ^ lO^g/cm^, whereas case b, based on l/V'- 
Sommerfeld-Chandrasekhar expansions ~ 0{^^^) holds for 
■0^5 (see Fig. 1 in paper I). 

Carrying out the expansions (~ 0{X'^)) in case a we get for 
Eqs. (0) and (||) 



Pc 



with the correponding particle number density 



h/2W ^"(V'.A) 



(28) 



(29) 



A3 



^i/2(^,A) + -AJ3/2(?A,A) 



(25) 



The thermodynamical potentials of the electrons are character- 
ized by the relativistic Fermi-Dirac integrals J,y{ip, A), which 
depend on the degeneracy {ip) and the relativity (A) parameters 



^ kT 



A 



kT 



(26) 



Generally, the degeneracy parameter ijj in ( |26| ) is a function of 
the density and temperature defined by Eq. (^s]). In order to 
evaluate the free energy 



(27) 



we have to determine ip explicitly from Eq. ( |25| ) by an inversion 
procedure which has to be performed numerically. Note, that tp 
in ( p6| ) is identical with the ideal contribution of g in (|l8|), i.e. 
V' = gl'^. Lamb (1974) and Lamb & Van Horn (1975) evaluated 
the thermodynamical potentials applying the parametrizations 
of Eggleton et al. (1973) for the relativistic Fermi-Dirac inte- 
grals. Johns et al. (1996) improved the accuracy of the poly- 
nomials given by Eggleton et al. (1973). Straniero (1988) cal- 
culated the complete set of thermodynamical functions based 



with the abbreviations L/^ and (not to be confused with the 
internal energy or the volume) 



[/^ = 1 
and 

= 1 



— A^' 



5/2 



'3/2 



1 



1a^ 

16 /. 



7/2 
5/2 



1 



3 /< 



-A- 



9/2 
I7/2 



(30) 



3/2 

/2 



7 /. 
16 I' 



5/2 
3/2 



3 /: 



7/2 

5/2 



(31) 



taking into account relativistic corrections. Ii, 



J^dz z'^(e^-'^ + l)-Vr(iy + l) are the nonrelativistic Fermi- 
Dirac integrals considered by parametrizations and expansions 
(see e.g. paper I). 

For the Helmholtz free energy in Eq. ( p7| ) we have per- 
formed the inversion ip — ip{nc, T) analytically (see paper I). 
Using the Maxwell relation Eq. ( pj[ ) and Eqs. (0)-(0) we obtain 
for the internal energy, compressibility, and coefficient of strain 



15 4/2 Ul 
8 h/2 V- 



1 _ h/2 V 



jLid 



1-1/2 



(32) 



(33) 
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3 I1/2 



2 1-1/2 



(34) 



Obviously, Eq. ( p2| ) results for A ^ in the well-known non- ya _ ^ 
relativistic relation Calculating the derivatives of 

the free energy with respect to the temperature in order to get 
the isochoric specific heat (Eq. [l^) we have to execute ex- 
plicitely 



tpX 



15^ I1/2 



O 1-1/2 

15/-1/2 
8 1-3/2 



7 I- 
1 + T7:A ■ 



7. /; 



1 + -A 



3/2 
/1/2 



3/2 
16 /1/2 



1 



5/2 
3/2 



8 /; 



5/2 



3/2, 



3,/: 



3/2 



8 h 



/2 



(43) 



(44) 



(45) 



-T 



(9r2 dT d'^dT 



+ 



d0 
dT 



9^ 



dr2 



.(35) 



The temperature derivatives of ?/; at constant V and A^e must be 
calculated from Eq. (Ey) giving 



dT 



3 /1/2 



2 /_ 



1/2 



4 4/2 ^4 



(36) 



From these relations we can read off immediately various 
asymptotics. The nonrelativistic (A <C 1) and nondegenerate 
limit (V" ^ —1) of Eq. results to the well-known classical 
value 

3 



(46) 



Furthermore in the nonrelativistic limit and at arbitrary degen- 
eracy the comparison of Eqs. (^^ and (38 ) yields 



(47) 



dT2 



j.2dV 



15 h/2 

4 1-1/2 



3 I1/2 /-3/2 



5 1-1/2 V- 1-1/2 V- 



I1/2 v. 



\1p 



h/2 1-3/2 



1-1/2 

h/2 



1-1/2 1-1/2 



h/2 ^aV 



-1/2 



-1/2 



I h/2 1-3/2 



V: 



ipip 



7 h/2 V, 



2/1/2 1-1/2 



30 h/2 



<37) 



Finally, Eq. (35) divided by N^k yields the ideal part of the 
electronic specific heat 



-V,c 



15 

T 



13/2 



3 h 



/2 



/1/2 V- 



ya 



5 /-1/2 V- 



105 , 



5/-1/2 V- 



h/2Ul_i h/2 V, 
h/2 V- 

h/2 

/1/2 V- 



15 /; 



3/2 



l^A h/2 IT 



7 /-1/2 h/2 V- 



where the abbreviations are given by 



7 /• 



[/J = 1 -f -A 



7/2 



1-5/2 



16 I' 



9/2 



7/2 



9^/9/2 

8 ^7/2 



7 L 



5/2 



3/2 



16 /. 



7/2 
5/2 



^8) 

(39) 
(40) 
(41) 



Case b describes the regime where the Sommerfeld- 
Chandrasekhar expansion (Chandrasekhar 1939, Cox & Giuli 
1968) becomes valid. The degeneracy parameter 1/; is related 
for the strong-degenerate limit by 



ip\ = Vi+^ 



1 



(48) 



where a = p-p/mc is the fraction of the Fermi- to the rel- 
ativistic momentum. The Sommerfeld-Chandrasekhar expan- 
sion with - 0{Tp-'^) provides (Cox & Giuli 1968) 



Pc 



4A V^{il;,X) 



(49) 



2 

A3 



2 a 
n3X 



with the abbreviations 



1 - 



3 

2^^ 

2^2A2 



1 - 



a" 

7 772 A2 



In {a + Vl + a^) 



60 



'l-2a^ 



(51) 



and 



7r2A2 
2a4 



1 



2a^ + 



7 TT^A^ 

20 a4 



(52) 



For the Helmholtz free energy represented by Eq. ( P7| ) the in- 
version V' = i^{nci T) for case b is already given by Eq. (p^. 
In order to correct Eq. (Q) to the order ~ 0(?A^^) analyti- 
cal inversions have been performed by Yakovlev & Shalybkov 
(1989) and paper!. 

For the internal energy, compressibility, and coefficient of 
strain we get 
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w 



1 

A 

1 
A 



(54) X„ = a—X 
oa 



(55) 



X^ 



a\——X 
oaoX 



(62) 



(63) 



Eq. (pOb yields by neglecting the order -if) ^ (Chandrasekhar 



Neglecting the order tl)^^ in Eq. (|3|) we get the Chandrasekhar IQSQJJhezer et al. 1986, Yakolev & Shalybkov 1989) 
approximation for the internal energy (Chandrasekhar 1939, 
Cox & GiuH 1968, Lamb 1974, Eliezer et al. 1986, Yakovlev 
& Shalybkov 1989) 



1 



mc fmc 



vy/l + ^n + 2a^) - — - In (a + Vl + c 

(VT+^il + 3a^) - (1 + 2a^))] (56) c^,, = ^ 
6 a \ / 2w 



which delivers the strong-relativistic (a :s> 1) 
and the nonrelativistic (a <C 1) 

2 



(64) 



(65) 



(66) 



The specific heat Eq. (|l5]) can be calculated for case b via 



Cv = -T 



d^F da d^F 

L O 

QT2 



dT dadT 
daV d^F 



limiting laws. Let us remark that both Eqs. ( |38| ) and ( |60| ) have 
in the degenerate, weak-relativistic limit a wide overlap, which 
we use to combine case a and b in order to realize numerical 
continuity (see paper I). 



dF 



dT da^ dr2 da 



-I V,Na 



^^"^^ 4.2. Exchange 



The temperature derivatives of a at constant V and must be 
calculated from Eq. ( pO| ) 



The lowest order exchange (Hartree-Fock) free energy for rel- 
ativistic electrons is given by 

c2 



T da _ 

adT~'3V^ + 



(58) 



J CO 



1 



(67) 



T2 



Tda 
a dT 



a dT 



T da 



a dT 



Finally, Eq. (p7h divided by Nek yields 



where J^i(V', A) is the relativistic Hartree-Fock integral 
(Kovetz et al. 1972, paper I). 

The A-expansion in case a results in 

J cc 



59) 



(68) 



where 



„Tda 
a di 



-'1/2 



1 V- 



° il/2 



AA 



a2 - 1 



yb 



Tda\ 

a dry 



6 (vTTa2 



1 1 



yb 
yb 



Ha^x 

72 /i/2 
A/ 



1 _ 25 ^^ /5/2 ^1/2 
1 1 -^3/2 -^3/2 



.(69) 



Parametrizations for the nonrelativistic Hartree-Fock integral 



1 

* act 
b" 



yb 4 yb 
T2 d^a 



(1 + ^2)3/2 



yb 



6 y 

yb 

_ 

yb 



(70) 



a dr2 



1 + a2 _ 1 



yb\i 



are given by Perrot & Dharma-wardana (1984) and Kalitkin & 
Ritus (1986) (see Stolzmann & Blocker 1996b). Note, that for 
A = the following relation holds 



Here the abbreviations labeled by the superscript b are given The exchange pressure pl^ = g^^ - 
by the simple derivatives with X ^ {U^,V^} exchange Gibbs energy (cf. paper I) 



^A = X-^X 



Q2 

^aa = A2_X 



(61) 9c 



kTK 



-1/2 



(71) 

/^(, is determined via the 
(72) 
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where 



1/2 



16 



1/2 



-1/2 



-^3/2 -^-1/2 
-^1/2 -^1/2 



24 / 



3/2 
1/2 



1 



25 /5/2 -^-1/2 



3 /; 



3/2 -^1/2 



(73) 



For the exchange contribution of the internal energy, compress- 
ibility, and coefficient of strain we get 



2Pc 



1 I- W-^ 15 /3/2 



2 kTK,h/2 



W, 



(74) 



(75) 




-'1/2 -'-3/2 
/-1/2 /-1/2 

,5, h/2V^ 



h/2 



15 /1^^ 

4 /-1/2 v; V- 



(76) 



Using ( |35| ) we obtain for the exchange part of the isochoric 
specific heat 





'15 






T 


M2 \y^~ 




3 

+ 5 


/1/2/-3/2 




^-1/2 



/3/2/-1/2 Vx 
/1/2 v" 

3 /3/2/-3/2 



'1/2 



V:,. 



225, 
"64"' 



^1/2^-3/2 T/;,^ 



/1/2/-1/2 y^i' 



I ^3/2^-3/2 



1/2 



/1/2 



^-1/2 ^ ^ 




with the relativistic corrections 
3 



/2 

w^f 1 + rA ^ 

A 8 jx 



1 



27/3/2 

8 /1/2 



30 /3/2 l^f 



216 h 



'1/2 



1 - 



25/ 



5/2 /1/2 



11 /3/2 /3/2 



(78) 



1 P 

1 ^1/2 



8 /-1/2 
9 ^^ -^3/2 



27 /; 



1- — A 



3/2 



'1/2 



1 + Ia'^ 



3/2 



72 / 



1/2 



1 - 



25/ 



5/2 /1/2 

11 /3/2 /3/2 



(79) 



1 + -^-1/2 -^-1/2 
/-3/2 

■1 + 3^ 

-'-3/2 -'3/2 



/1/2 

/1/2 7 /-1/2 

— A — X 

9 /_3/2 



16 /1/2 

(l -|_ 3 /1/2 /1/2 25 4/2 /-3/2 " 
V 28 /_i/2 /3/2 28 /3/2 /-1/2 / 



(80) 



W: 



1 + -A 

4 /_ 



1/2 



1/2 



^27/3^ 
16 /1/2 

125^/5/2 



/1/2 /1/2 

/-1/2 /3/2 



216 /. 



3/2 



1 



3 / 



3/2 



25/ 



5/2 



(81) 



The Sommerfeld-Chandrasekhar expansion for the exchange 
free energy in case b (Salpeter & Zapolsky 1967, Kovetz et 
al. 1972, Lamb 1974) deUvers for - 0{iIj-^) (see paper I) 



•'^^ kTK 2a3 



1 W^\a,\) 



(82) 



where 



2(1 + a2) 
^2 



3 

TT A 
18 



-0.7046- 



2 In 



2a' 



3B 



10^2 



37 



a 

6W\ 



lla2 llofS / 



(83) 



with B = \/m?hi 



VT 



With the Gibbs energy (see paper I) 
e2 3 1 

9cc 



fcTAc 2a3 V2^ + 



(84) 

J ce 



it is possible to provide the exchange pressure p^^ — g 
again. 

In a straightforward manner, the exchange contribution for 
internal energy, compressibility, and coefficient of strain are 
represented by 



1 



kTA, 2a3 J2^ 



1/b 



( 



y^ 



y}^ 



(85) 



'''T,cc 



fcTAe 2a3 V3^'' + KS' 



yb 3^b ^ m, 



yb 



lyb - 

6 



(86) 
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"S,ee 



Pee + < 



9e. 



1 



The exchange specific heat can be formulated by 



(87) 



1 



a AT 



Tda 
a dT 



yb 



+ 



T2 d^a M^, 



a dT2 



8) 



The temperature derivatives of a in Eq. ( |88[ ) are given by i 
and (|9|). The derivatives of W with respect to A and a for the 
exchange contributions are determined by 



Wl = a—W"" 
oa 



^b ^ aX— — 
dadX 



AA 



9A2 



(89) 
(90) 
(91) 



The adiabatic temperature gradient (||) expressed by our dimen- 
sionless potentials ( [l8| ) and ( [l9| ) is given by 



Vad = Pfcx = —(1 ) 

Cp 0s cp 

with the isobaric specific heat 



Cp = cv 



(92) 



(93) 



The Figs. |l]and ^ illustrate the adiabatic temperature gradient 
(92) of the ideal gas im comparison with a gas where addition- 
ally the exchange term is taken into account. Note, that the rep- 
resentation in Figs. |l] and |^ for extremely high densities serves 
only to illustrate the asymptotics of the theoretical expressions. 
Processes relevant to such high densities as, e.g., pycnonuclear 
reactions or electron captures, are not taken into account. Fig. |] 
illustrates the course of Vad for the elements hydrogen, helium 
and oxygen along a nonrelativistic and a weak-relativistic tem- 
perature. A comparison with the corresponding data for carbon 
from Lamb (1974) is given in Fig. |[ Note a disagreement for 
T = lO^K. This deviation could be caused by an erroneous 
constant for the exchange contribution used by Kovetz et al. 
(1972) and Lamb (1974) and is discussed in detail in paper I. A 
numerical procedure to guarantee a smooth transition between 
the approximations given by case a and case b is developed 
in paper I. The curves as seen in Figs. [l| and ^ show how ex- 
change effects vanish with increasing temperature. Considering 
the isotherms in the high density region exchange effects are 
negligible. Nevertheless Eqs. (|8^)-(|8^) guarantee the accurate 
asymptotics, which is a fundamental problem for the calcula- 
tion of EOS derivatives for quantities with vanishingly small 
temperature dependence (Miralles & Van Riper 1996). 




4 6 8 10 
log [Q/{g cm"^)] 

Fig. 1. Adiabatic temperature gradient Vad vs. density at T = 
lO^K and T — lO^K for various elements. The curves refer to 
the sum of ideality, exchange and radiation. The presentation 
towards extremely high densities serves only to illustrate the 
asymptotic behaviour of the adiabatic gradient in our model 
(processes relevant to very high densities as, e.g., pycnonuclear 
reactions or electron captures, are not considered). 

4.3. Correlations 

According to Eq. ( |20| ) we have to determine the terms labeled 
by the superscript c, which are based on correlation effects. 



0.50 - 




2 4 6 
log [Q/{g cm"^)] 

Fig. 2. Adiabatic temperature gradient Vad vs. density for car- 
bon along various isotherms compared with Lamb (1974) (+ : 
lO^K, A : lO^K, • : lO^K, o : lO^K, * : lO^K). The curves 
refer to the sum of ideality, exchange and radiation. The presen- 
tation towards extremely high densities serves only to illustrate 
the asymptotic behaviour of the adiabatic gradient in our model 
(processes relevant to very high densities as, e.g., pycnonuclear 
reactions or electron captures, are not considered). 
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As in paper I we aim to present closed-form parametrizations 
formed by Pade Approximants to get explicite expressions for 
the correlation contributions of the thermodynamical poten- 
tials. The thermodynamical potentials resulting from the cor- 
relation contributions are characterized by the electronic and 
ionic Coulomb-coupling parameters To = (47rne/3)^/^e^/fcT 
and Fi = TdZ^^"^) giving the strength of the Coulomb in- 
teraction, and by the degeneracy parameter 8 = T/Tp = 
2(4/97r)^/'^rs/ro describing the quantum state of the system. 
The parameter rg = nioe^ /h^ (S/ATmc)^^^ is the ratio of the 
mean interelectronic distance to the (electronic) Bohr-radius. 

4.3.1. Electron-electron interaction 

The Pade Approximants for free Helmholtz- and Gibbs ener- 
gies of the nonrelativistic electronic subsystem at arbitrary de- 
generation and Coulomb coupling have the general structure 



J cc 



9c 



3/2 



1 + aiVl^^ + 
S2r^ [Mc(rs,0) + AMc(rs,r)]/T 



1 + Sir; 



3/2 



S2 



(94) 



(95) 



with ground-state energies ec{rs,0) and /ic(rs,0) derived by 
Vosko et al. (1980) and low-temperature corrections Aedrs, t) 
and A/ic(rs,T), which are given in paper I. The coefficients 
read 



o-a — —^foiTc) , 



/o(r) 



3 1 dao 

Sn = —an H i o 

" 2 " 3 dTc 



1 



ai 



a2 



3V3 
"32"' 



(1 + 0.1088re)i/2 (1 + 0.3566re)'V2 
277 



'2'KT 



S2 = 0-2 



(96) 
(97) 

(98) 
(99) 



Note that we have modified the coefficients oi and si as com- 
pared to paper I by a rearrangement of the quantum virial func- 
tion (EbeUng et al. 1976) 



Kl(x)=E*{-x) 



8x 



Ql{-x) 



ln|a;|-Co 



In paper I we incorporated for (x) in Eq. ( |98[ ) only the first 
part of Eq. (100) as proposed by Ebeling (1993). The second 
term was considered by the coefficient ci of the ion-electron 
contribution (see Eq. (85) in paper I). The justification to rear- 
range these summations over the electrons is the omission of 
the coefficient as in Eq. (53) in paper I, which has only been 
introduced to optimize the interpolation. For E2 and Q3 we 
derived Pade approximations as given in paper I with improve- 
ments to achieve higher accuracy. 



El{-x) 



ln2 + 0.0113a;2 + 0.1a;' 



l + 7r3/22;/(181n2) + 0.1x5 




Fig. 3. The isochoric specific heat for the electronic subsystem 
at different degeneracy (0 = 10 (classical), 1 (intermediate), 
0.1 (strong)) vs. Coulomb-coupling parameter. 



Qli-x) = 
with Co 



(\\n\x\+Ca 



8a; 



1 



lOCi 



\^ 192 



75 



(102) 



i(ln3 + 2CE-^) and Ci 



6 V " ' 6 

i(ln3+%-i). 

A detailed comparison of the Pade formula (M) with most 
advanced many body calculations is given in paper I and re- 
cently in Stolzmann & Rosier (1998b). The electronic correla- 
tion contribution for the pressure is obtained by 

The internal energy can be described by the Pade Approximant 



u^vl'^ - U2TI [ec(rs, 0) + Aee(rs, r)] /r 



3/2 



(100) ui 



with the coefficients 

, doo 
"dTe 

r dfli 
3'd7 



U2 = 02 



(104) 



(105) 



^_a^~^^(r^K*{x)]{m6) 



The temperature derivative of the quantum virial function (|10 
is given by 



r-Kix) 



r^E*{-r)-\K:{-r) + \E;{-r) 



16 



/2tit 



d ^ 1 



dr 



12 



(107) 



For the temperature derivatives of E2 and Qg we derive the 
Pade Approximants 



(101) ■r^E;{-T) = 



36 12 \^ 



dT 1 _l_0. 113 J%a;(l- 0.473a;) + f| 



(108) 
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dr 



384 - 



O.lx 



4.5 



16a; 



1 



384 , 

2571^2' 



0.1a; 



4.5 



.(109) 



The correlation contribution l/fcx co fo'" '^^e bulk modulus is 
given by 



3/2 



fc2r6[pe(r„0)(l-l/fcc(rs,0))]/7 



^T,ee 



1 + kiV. 



3/2 



with the coefficients 

3 Tc dao 

^o = -ro-Tdf: 



fc2 = 10a2 



The ground-state bulk modulus is given by 

1 = ? !f A„ ir 0) 

fcc(rs,0) 0.062181 3 dr/"^ " ^ 



(110) 



/ci = 30ai .(111) 



(112) 



4.3.2. Ion-ion interaction 

The classical one-component plasma (OCP) is a thoroughly 
studied plasma system and the results are unified in closed- 
form parametrizations (Hansen 1973, Graboske et al. 1975, 
Brami et al. 1979, EbeHng & Richert 1985, Ichimaru 1993, 
Kahlbaum 1996, Chabrier & Potekhin 1998). For the ion-ion 
correlation we modify the Fade Approximants used in paper I 
by the new formula (Stolzmann & Blocker 1998, 1999) for the 
free energy 



•I 11 



3/2 



1 + bzVl'^F{Ti) 



62rfeii(ri) 



1 - 6irf G(ri) + b^Tf 

with F(ri) = InTi + Bq and G{Ti) = InTi + Bi 
3 



So 



Using the Monte-Carlo data-fit from Vosko et al. (1980) for 
Pelt's, 0) (see paper I) we apply 



2Ce 



■ln3 



1 



fcc(rs,0) 



3/2 



l + miy^ 2 + TOiy^ + m2r 



3/2 



Bi = 3 I 2Ce 



iln3- 



11 



21n2- 



17 

— I - 0.4765 



(113) 

and the coefficients 



The coefficient of strain 01 can be expressed by 



6o = 



(Z2)3/2 



P. 



V3(Z)l/2(Z5/3)3/2 



hi 



3V3 



bo = 100 



-.3/2 



9li c 



with the coefficients 
9 



16 



-ai 



3 dap 

4'Mre 

2 dfli 
' 3^ d7 



2«2, 



(120) 

(121) 
(122) 

(123) 
(124) 



Eq. (12C) is based on the classical Fj < 1 result from Cohen & 
Murphy (1969) and for Fj ^ 1 we take into account the most 
(115) recent Monte-Carlo fit for the free energy £ii(Fi) of the liquid 
OCP from DeWitt & Slattery (1999) 



(116) = 0.899172Fi + 0.274823 InFi 

-1.864179F;°-323064 



1.4018 .(125) 



and for the isochoric specific heat Cy the Fade Approximant 
can be obtained by 



We remark, that in our notation of Eq. (g^ for £ii(Fi) we have 
to take for Fi ^178 



wqf; 



,3/2 



with the coefficients 

3 " 



[2/(r,)e2(l + |lne)] /r 



1 + t^iF 



3/2 



V2 



Vl 



ao 



ir 

4'^dFe 



V2 = 02 



(117) eii(Fi) = 0.895929Fi- 1.51nFi + 3.9437Fri 

+ 1245Fr2 + 1.1703 

which, with Eq. ( pT] ) rewritten by 

(118) o 

/i'=ol^^ + 31nFi- 0.7155, 
2 Kydi 



(126) 



(127) 



3%/3 
32 ' 



'2ttt 



+ ^(l + ^2(-r)) 



5 16 



6 x/2^ 



V2t:t 
16 

d 

d7 
1 

12 



Ki-r) 



delivers the thermal (solid phase) energy of the OCP given 
by Hansen (1972), Pollock & Hansen (1973), Slattery et al. 
(1980), and Stringfellow et al. (1990). The quantity Rydi = 
AZ'^mp/rricRyd denotes the ionic Rydberg energy with the 
(electronic) Rydberg unit Ryd = mae^ /2Tt? . 

The free Gibbs energy gf; and the isochoric specific heat 
^^^^^ Cy 2st cxpressed by the Pade Approximants 



Fig. ^ shows that the electronic correlation part of the isochoric 
specific heat can be neglected in the case of strongly degenerate 
electrons. 



91 



3 1. p3/2 
2°0l i 


l+4^3rf/2 (F(FO + i) 


+ fo2FfA*ii(Fi) 


l-|6iFf (G(FO + ^) 


;i + frp ; 


+ 62Ff 



(128) 
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10' 
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Si 



Fig. 4. The compressibility for the ionic subsystem vs. coupling 
parameter Fj. The solid line represents the correlation contri- 



Fig. 5. The free energy contribution of the ionic quantum cor- 
rections for carbon at T = lO^K vs. ionic quantum param- 



bution (133) and the dashed line refers to the parametrization 



from Hansen (1973). 



3l p3/2 



with 



1 + 863rf/' (i^(ri) + 1) + 662rf fti(ri) 



76irf (G(ri) 



32^ 
63 y 



-(12' 



eter (Eq. 135). The solid line represents Nagara et al. (1987) 
(Eq. 137) (shown only within its validity regime), the dashed 
line refers to Chabrier et al. (1992), and the dashed-dotted line 
to lyetomi et al. (1993). The curve labeled by WK refers to 

e first order in the Wigner-Kirkwood expansion and HV to 

ansen & Vieillefosse (1975). 



d^gii(ri) 



temperatures. A measure for the quantum effects of the ions is 
.(130)the parameter (weight fraction Xi, molecular weight Ai, proton 



mass TTip, atomic mass unit TOu = 1.66053 • 10 g) 



The pressure, internal energy, isothermal compressibility, and 
coefficient of strain (Hansen 1973, Stolzmann & Blocker 1999) 
are calculated by 




Pn 



U:; 



1 1 , 



^T,ii 



9 



^V,ii ■ 



Cv,ii 



(131) 
(132) 

(133) 
(134) 



which is related by 




(135) 



(136) 



Fig. ^compares our Pade Approximants of the isothermal com- 
pressibility with the fit formulae from Hansen (1973). Further 
comparisons between different interpolation formulae for the 
classically ionic subsystem, e.g. from Hansen (1973), Ichimaru 
(1994), Kahlbaum (1996), Chabrier & Potekhin (1998) can be 
found in Stolzmann & Blocker (1999). 

Ionic quantum effects, which were included approximately 
in our earlier Pade Approximants for the ionic subsystem (see 
paper I) will be considered separately in the following section. 

4.3.3. Ionic quantum effects 

The classical description for the subsystem of ionic OCP be- 
comes inaccurate in the region of high densities at moderate 



and whereas ojp denotes the (ionic) plasma frequency and 
Rs = TOiZ2e^/?i^(3/47rni)i/^ is the ratio of the (ionic) 
Wigner-Seitz -radius to the (ionic) Bohr-radius. Quantum cor- 
rections for a solid OCP over a broad region of Qi were calcu- 
lated first by Pollock & Hansen (1973). Chabrier et al. (1992) 
(see also Chabrier 1993) and lyetomi et al. (1993) derived an 
free energy for a solid OCP model with quantum coiTections. 
Nagara et al. (1987) calculated quantum corrections for the free 
energy of the OCP valid for fluid and solid phases by a new ex- 
pansion method. In our model ionic quantum effects according 
to Nagara et al. (1987) as well as Chabrier et al. (1992) wifl be 
included. 

Following paper I, the free energy calculated within the so- 
called "sixth reduced moment approximation" of Nagara et al. 
(1987) is given by 



3 t 

qo{t) + Ziqi{t) - - — 2595(0 + zeqe{t) 



(137) 
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Fig. 6. The free energy contribution of the ionic quantum cor- 
rections at Fi = 180 vs. 9i. The solid line represents Nagara 
et al. (1987) (Eq. 137) (shown only within its validity regime), 
the dashed line refers to Chabrier et al. (1992), and the dashed- 
dotted line to lyetomi et al. (1993). 

with t ~ 8i/2\/3. Generally, the reduced moments Zi are de- 
pendent of Fi. For simplicity we choose for the solid phase the 
values for the rigid bcc lattice (Nagara et al. 1987): 



Zi = 0.827702 , Z5 = 1.131 
The functions qi (t) are given by 
/ sinht^ 



ze = 0.55045 



qo{t) = In 



94(i) = ^ - -cothi 



1 



sinh t 



/ N ^ , 2 3 , 5 1 
95 (i) = o coth t - - cothi - -t n 
■3 Z b Sinn t 



2 

+ 7 



(138) 

(139) 

(140) 
(141) 



1 * , 
D lb 



1 



16 sinh-" t 



1 



2t 



cothi .(142) 



The free energy expansion ( |137| ) reproduces the Wigner- 
Kirkwood expansion (Wigner 1932, Kirkwood 1933) up to the 
order 8f . Furthermore the expansion of Nagara et al. (1987) 
delivers the high-0i asymptotics of the Coulomb lattice by 



(143) 



with fii — 0.564. For the harmonic lattice approximation one 
obtains ni = 0.511 (Hansen & Vieillefosse 1975). 
The ionic quantum effects for our set of potentials are given by 
the resulting formulae 



cq 
Pii 



-t 



3 t 



3 t 



- 



1 1 


1 1 1 1 M 1 III 


II 1 1 1 1 1 1 1 






log p = 7.5 


- 
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Fig. 7. The isochoric specific heat for the ionic quantum cor- 
rections for carbon at p = 10^'^g/cm^ vs. temperature. The 
solid line represents Nagara et al. (1987) (Eq. 147) (shown only 
within its validity regime), the dashed line refers to Chabrier 
(1993), and the crosses to Pollock & Hansen (1973). The curve 
labeled WK refers to the first order in the Wigner-Kirkwood 
expansion and HV to Hansen & Vieillefosse (1975). 



^ 2 

\ 



T 



log Q-7 .5 
Carbon 




10' 



10' 



T/K 



Fig. 8. The electronic and ionic parts of the isochoric specific 
heat for carbon at p ^ lO^'^g/cm"^ vs. temperature. The long- 
dashed line represents the electronic terms of Eq. ( |20| ) The solid 
line (applies Nagara et al. (1987) for c^jj) and the short-dashed 
line (applies Chabrier (1993) for c^jj) represent the ionic terms 
of Eq. (|y). The crosses refer to the ionic contribution calcu- 
lated by Lamb & Van Horn (1975). 



cq cq 
9n =Pii 



f- 
•I 11 



cq 



3 t 



(144) y::^^2^^ + --z,q,{t) 



(145) 



(146) 
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-V,ii 



l,cq 
'*'T,ii 



3 t 
if; 



cq 



(147) ''V.ii 



S,ii 



(158) 



cq 



3U( 



-V.ii 



^V,ii 



3 t 

25 

4ri 



(l5{t)-tq'S) 



(148) 



(149) 



The functions q[{t) and are first and second derivatives 
of the functions qi [t) with respect to t. 

Although the expansion of Nagara et al. (1987) can be ap- 
plied to estimate quantum effects for fluid ionic phases as well 
as for solid phases this expansion is restricted to <di ^ 17, 
which is equivalent to Ti 4 • 10^{AZ'^/{T/K))'^/^ or to 
p <, 5 ■ 10~^(T/K))2(A/Z)2 for the density-temperature 
regime. For Fi w Fm <; 180 it is convenient to use the quantum 
crystal models from lyetomi et al. (1993) or from Chabrier et 
al. (1992), which are applicable for a wide region of Qi. 
If we extract the ionic quantum corrections for the free energy 
of a crystal from Chabrier et al. (1992) one obtains (Stolzmann 
& Blocker 1994) 



fi^ = -3i?3(«eO + 21n[l 



In [l 



+/3ei-31n {6Qi) 



with (Chabrier 1993) 



a = 0.399 



7 = 0.899 , (3 = 0.767 



(151) 



The last term in (150) with S = 0.4355 guarantees the pure 



ionic quantum correction in our formalism. (77) is the De- 
bye integral, which can be approximated by (see e.g. Iben & 
Tutukov 1984) 



V Jo - 1 



1 + 0.43?? - 



With (150) our set of potentials is given by 

C( 

Pii 



1 7eiexp(-79i) 3 1 

2 l-exp(-7ei) 2 2^ ' 



an"" 



cq 



rcq 
J ii 



2Pu 



cq 



.(152) 



(153) 



(154) 



(155) 



An explicite expression for (158|) is given by Chabrier (1993) 
Figs 



compare the ionic quantum corrections for the free 
energy as well as the isochoric specific heat calculated by dif- 
ferent approaches. Note, the end of the solid lines in the Figs. ^ 
P reflects the validity of the moment expansion as pointed out 
by Nagara et al. (1987). Fig. | shows ionic quantum correc- 
tions for Z = 6 calculated with different theories. Ionic quan- 
tum corrections in the vicinity of phase transition at Fj « 180 
are shown in Fig. ^. Fig. ^ illustrates that the quantum cor- 
rections calculated by lyetomi et al. (1993) decrease with in- 
creasing Gj. This behaviour is caused by the inclusion of the 
anharmonic contribution of the zero-temperature oscillations 
(see e.g. Nagara et al. 1987, lyetomi et al. 1993), which is ne- 
glected by Chabrier et al. (1992). Fig. ^ shows the ionic quan- 
tum corrections for the isochoric specific heat derived from the 
corresponding free energies. Fig. ^illustrates the course of the 
electronic and ionic contributions of the isochoric specific heat 
vs. temperature and shows how different prescriptions for the 
quantum corrections affect the results. The discontinuities in 
the curves for the ionic contribution in Fig. ^ refer to the fluid- 
solid phase transition at Fi « 180. 



(150) 4.3.4. Ion-electron interaction 



Next we have to include the screening effect between electrons 
and ions. We adopt the new Fade Approximant based on the 
Unsold-Berlin-Montroll asymptotics of Stolzmann & Ebeling 
(1998) and EbeUng et al. (1999) 



coFf^^ + C2Fo'^eic(rs, Fi) + C3C4F; 



-.5/2 



f.c 

J IC 



5/2 



1 + ciVf^ + C2F; 



-^5/2 



C3F; 



5/2 



(159) 



The weak coupling limit (Fj <C 1) is given here by the Debye- 
Hiickel two component plasma (DHTCP) law. The coefficients 
Ck are defined by 



Co 



Cl 



1 (Ce + (Z2))3/2-Ce'/2_(Z2)3/2 



V3 



2co(Z)(Z5/3)3 



CaCbZlZl 

ah 



(160) 



8 



Q^{-\U\) 1, I, 



I -Co 



(161) 



^T,ii 



^1% 



3 a6iexp(— a0i) 1 
2 1 - exp(-aei) ^ 4 



"-T.ii 



1 76iexp(-79i) 
4 1 - exp(-76i) 
3 a6iexp(— a6i) 
"2 1 - exp(-Q:ei) 

, 1 



76; 



1 - exp(-7ei) 



(156) 
(157) 



C2 

C4 



Zr-'^l^ , C3 = 0.01 , 

(2- 1.786- 0.9 - 0.9Fe/Fi)ri 

2-Za-^h 



VT(7a + 7b) ' 



Ik 



mo . 
nik 



(162) 
(163) 

-(164) 



In Eq. ( |161[ ) is summed up over a 7^ 6 in contrast to sum- 
mations made in earlier approaches (Ebeling 1993, paper I). 
For eio(r's,Fi) we apply the expression given by Ebeling & 
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Fig. 9. The free energy contribution of the screening con- 
tribution for hydrogen and densities referring to — 
2,1,0.5,0.1,0.01 (from top to bottom) compared with the 
Debye-Hueckel li mitin g law (DHTCP) (the first term in the 
nominator of Eq. (159)) and the high density /large coupling- 
asymptotics from Galam & Hansen (1976) (dashed-dotted 
lines). The dashed lines refer to Chabrier & Potekhin (1998) 



Fig. 11. The free gibbs energy contribution of the screen- 
ing contribution for carbon and densities referring to = 
2,1,0.5,0.1,0.01 (from top to bottom) compared with the 
Debye-Hueckel limiting law (DHTCP) and the high den- 
sity/large coupling-asymptotics from Galam & Hansen (197 6) 
(dashed-dotted lines). The solid lines correspond to Eq. (167). 



and the solid lines correspond to Eq. (159) 



asymptotics from Galam & Hansen (1976) including high- 
density corrections according to Kagan et al. (1977). 



Richert( 1985), which consider the high density/large coupling- eie(rs,ri)= j Tir^ [e['c(?'s, Ti) + J(r)] (165) 



12Z\ 



2/3 



B r: 



-3/4 



0.5 



0.0 



-0.5 



o 



-1 .5 



-2.0 



-2.5 




A + Crr 



1 + 0.089r2 1 + 0.353 r: 



-i/i 



(166) 



with the relativity parameter x = 0.014/rs. J(ri) as well as 
A, B and C depending on Z are given by Galam & Hansen 
(1976). Recently, Chabrier & Potekhin (1998) presented a 
Pade formula for the screening contribution of the free energy, 
which includes the Thomas-Fermi approximation in the limit 
of high densities and large coupling parameters. The results of 
Chabrier & Potekhin (1998) with Eq. (159) are compared for 
hydrogen and carbon at five densities in Fig. ^ and Fig. |o[ 
In the following we proceed in the same manner as in the case 
of the electron-electron subsystem in Sect. [4.3. 1[ 
For the ion-electron Gibbs energy we use the Pade expression 



5i< 



C2ro^^ (eic + \ {hie - die)) + C3C4T. 



5/2 



1 



Cirf/'+CaFe'^+Cgn 

with hie = Ti • deic/dVi and die = ■ dsic/drs. Fig. [TI| 
illustrates the Gibbbs energy for carbon at five densities. The 
polarization pressure can be determined by p?, = gf^ — /j^. 

5. Numerical results and discussion 

This section will be devoted to illustrate the behaviour of the 
adiabatic temperature gradient Vad- The adiabatic temperature 



5/2 



-.5/2 



-(167) 



Fig. 10. The free energy contribution of the screening 
contribution for carbon and densities referring to rg = 
2,1,0.5,0.1,0.01 (from top to bottom) compared with the 
Debye-Hueckel limiting law (DHTCP) and the high den- 
sity/large coupling-asymptotics from Galam & Hansen (1976) 
(dashed-dotted lines). The dashed lines refers to Chabrier & 



Potekhin (1998) and the solid lines correspond to Eq. (159) 
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Fig. 12. Contributions of the isochoric specific heat in units of 
krii for helium at T = lO^K vs. density. The lines refer to 
the ideal electrons (dashed dotted) and ions (long-dashed), the 
exchange and correlation (dashed), the ionic quantum correc- 
tion from Chabrier et al. (1992) (short-dashed), and the sum 
of those (solid). The curves with dots represent ionic quantum 
corrections after Nagara et al. (1987). 



gradient is a highly sensitive quantity, because it depends on 
first and second-order derivations of the model Helmholtz free 
energy. Therefore it is well suited to test the quality of the ex- 
plicit expressions derived in Sects. 3 and 4. The adiabatic gra- 
dient is a crucial ingredient for many applications. For instance, 
Vad determines the Schwarzschild criterion for convective in- 
stability and an accurate evaluation is required for evolutionary 
calculations of very-low mass stars. We restrict our compar- 
isons of the adiabatic temperature gradient to the EOS tables 
given by Straniero (1988), Saumon et al. (1995), and Rogers 
et al. (1996). Further detailed comparisons can be found in 
Saumon et al. (1995) who compared their results for Vad with 
those of Fontaine et al. (1977), Dappen et al. (1988), and Magni 
&Mazziteni(1979). 



Fig. 12 shows the total isochoric specific heat according to 
Eq. ( pO| ) as well as their various single contributions (in units 
of njfc), i.e. electronic exchange and correlation terms (elec- 
tronic, ionic and screening part). We note that the discontinu- 
ity in the ionic correlation is caused by the fluid-solid phase 
transition at F = 178 (Stringfellow et al. 1990), which can 
also be observed for the adiabatic temperature gradient in the 
high-density region. The influence of quantum effects on the 
fluid-solid transition of the OCP is discussed in detail by Pol- 
lock & Hansen (1973), Lamb & Van Horn (1975), Nagara et al. 
(1987), Chabrier (1993) and lyetomi et al. (1993). We chose for 
our numerical calculations F = 178, which is based on an sim- 
ulated OCP without quantum effects (Stringfellow et al. 1990). 
The specific heat contribution according to the ionic correla- 
tion is monotonously increasing until the phase transition of the 
OCP model. The first maximum at logp « 1.5 in the dashed 



log tp/Cg cm"^)] 

Fig. 13. High-density asymptotics of the adiabatic temperature 
gradient Vad vs. density for hydrogen, demonstrated along the 
isotherm T ~ 5 ■ 10"^ K. The dashed line refers to the ideal and 
exchange terms only. The dashed-dotted (without ionic quan- 
tum effects), the solid (with ionic quantum corrections from 
Chabrier et al. 1992), and the short-dashed (with ionic quan- 
tum corrections from Nagara et al. 1987) lines take into account 
correlations. For more explanation: see text. 
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Fig. 14. Adiabatic temperature gradient Vad (without radiation 
and all terms in Eq. (|2C|)) for hydrogen (solid line) at T = 
j^q5.54j^ vs. density. The dashed line refers to Saumon et al. 
(1995) and the dashed-dotted line to Fontaine et al. (1977). 



curve in Fig. [T^ is determined essentially by the electronic ex- 
change and electronic correlation. The influence of electronic 
exchange and correlation increases with decreasing tempera- 
tures. Quantum effects for the specific heat according to the 
ions after Chabrier (1993) are illustrated by the dotted line in 
Fig. |l^, which can be seen also in Figs. ^ and ^. This contribu- 
tion is compensated at very high densities by the contributions 
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Fig. 15. Adiabatic temperature gradient Vad (without radia- 
tion and all terms in Eq. (^) for helium at T = 10^ ''*K 
vs. density. For the ionic quantum corrections the results from 
Chabrier et al. (1992) and from Nagara et al. (1987) (soHd Hne) 
are used. The dotted line refers to Saumon et al. (1995) and 
the dashed-dotted line to Fontaine et al. (1977). The ideality 
(long-dashed line) is shown for comparison. 
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Fig. 17. Adiabatic temperature gradient Vad (with radiation 
and all terms in Eq. (pO|)) for carbon (solid line and short- 
dashed lines) at T = 10^ vs. density. For the ionic quantum 
corrections the results from Chabrier et al. (1992) (solid line) 
and from Nagara et al. (1987) (short-dashed line) are used. The 
triangles refers to Straniero (1988), the dashed-dotted line to 
Fontaine et al. (1977), and the crosses to Lamb (1974). 
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Fig. 16. Adiabatic temperature gradient Vad (with radiation 
and all terms in Eq. (^0|)) for helium at T = lO^K vs. density. 
For the ionic quantum corrections the results from Chabrier et 
al. (1992) (soHd Une) and from Nagara et al. (1987) (short- 
dashed line) are used. The dotted line refers to Saumon et al. 
(1995), the triangles to Straniero (1988), the dashed-dotted Hne 
to Fontaine et al. (1977), and the dots to Rogers et al. (1996). 
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Fig. 18. Adiabatic temperature gradient Vad (with radiation 
and all terms in Eq. (pO|)) for carbon (solid line and short- 
dashed lines) at T = lO^K vs. density. For the ionic quantum 
corrections the results from Chabrier et al. (1992) (solid line) 
and from Nagara et al. (1987) (short-dashed line) are used. The 
triangles refers to Straniero (1988), the dashed-dotted line to 
Fontaine et al. (1977), and the crosses to Lamb (1974). 



of the ideal and the correlated ions. Considering ionic quantum 
corrections calculated by Nagara et al. (1987) the isochoric spe- 



cific heat will be more strongly reduced as shown in Fig. 12 



We point out that the same behaviour appears for the adia- 
batic temperature gradient as indicated in Figs, [l^ and |l5[|l^. 



Note, that although the ionic quantum correction calculated by 
Chabrier et al. (1992) and Chabrier (1993) are valid only for 
the solid phase we considered here also the fluid range in order 
to illustrate the differences to the result of Nagara et al. (1987). 
Figs. [l3[]20| display isotherms of the adiabatic temperature gra- 
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Fig. 19. Adiabatic temperature gradient Vad (with radiation 
and all terms in Eq. (pO[)) for helium along the isotherms 
T = 10^ (solid), 10^ (dashed), and 10**K (dashed-dotted) vs. 
density. 
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Fig. 20. Adiabatic temperature gradient Vad (with radiation 
and all terms in Eq. (pO[)) for a mixture of C/O=50/50 (mass 
fractions) along the isotherms T — 10^ (solid), 10^'^ (dashed), 
and lO^K (dashed-dotted) vs. density. 



dient evaluated on the basis our analytic expressions given in 
Sect. Q including the photonic contribution (Figs. p"6[p0| ) over a 
broad range of densities. As already mentioned, the representa- 
tions towards very high densities aimed at the purpose to show 
the asymptotic behaviour of the potentials given by Eq. (pO[). 
The accurate description of plasmas at very high density re- 
gion requires the consideration of additional effects as electron 
captures and other processes (cf. Shapiro et al. 1983). Fig. 
shows the high density asymptotics of the adiabatic tempera- 
ture gradient calculated by different approximations. Consid- 
ering ideal and exchange contributions we obtained the well 



known value Vad = 1/2, according to relativistic consider- 
ations. The limiting value Vad — 3/8 is determined by in- 
clusion of the classical ionic correlation and Vad = 1/4 is 
obtained with ionic quantum corrections given by Chabrier et 
al. (1992). In Fig. p^p^ we compare our formalism with data 
from other authors (Lamb 1974, Fontaine et al. 1977, Straniero 
1988, Saumon et al. 1995, Rogers et al. 1996) based on analyt- 
ical as well as numerical methods and valid for various density 
and temperature ranges. The radiation contribution is dominant 
at high temperatures and low densities. The well-known limit- 
ing case of Vad — 0.25 in the low-density limit depends on the 
temperature in a strong manner. 

The largest deviations for the adiabatic temperature gradi- 
ent are observed in the intermediate density region, that means 
for plasma parameter where e.g. correlation effects must be de- 
scribed by accurate expressions. It should be emphasized that 
the Pade technique is well suited to meet the challenging re- 
quirements of this - sometimes poorly known - density regime. 

Fig. |l9| shows the adiabatic temperature gradient for he- 
lium. In Fig. ^we draw Vad for a mixture of carbon and oxy- 
gen at high temperatures as appropriate for white dwarf interi- 
ors. 

The temperatures and densities are chosen here according 
to the conditions for a fully ionized plasma. Nevertheless, our 
EOS-formalism can be applied for the conditions of a partially 
ionized plasma (lower temperatures) provided additional terms 



in Eq. (20) are considered and are solved in connection with 



the formalism of the dissoziation- and ionization equilibrium 
(Beule et al. 1999). 

6. Summary 

Our aim was to provide data for the thermodynamical poten- 
tials summarized in Eq. (|o|) which are applicable to, e.g., stel- 
lar interiors (fully ionized regions). For the Helmholtz free en- 
ergy, the Gibbs free energy, and the pressure (equation of state), 
expressions for each term in Eq. (pO|), i.e. ideality, exchange 
and correlation, have already been presented in paper I. In this 
paper we have improved our Pade Approximants for the corre- 
lation contributions and continued to present formulae for the 
explicit calculation of the internal energy, reciprocal compress- 
ibility (bulk modulus), coefficient of strain, and the specific 
heat. Further results using our calculational concept for the adi- 
abatic temperature gradient for fully ionized plasmas are pre- 
sented by Stolzmann & Blocker (1998). The extension of our 
formalism to the partially ionized region is in progress. 
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